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Abstract 

When rapidity gaps in high-p^ dijet events are identified by energy flow 
in the central region, they may be calculated from factorized cross sections in 
perturbative QCD, up to corrections that behave as inverse powers of the cen- 
tral region energy. Although power-suppressed corrections may be important, 
a perturbative calculation of dijet rapidity gaps in pp scattering successfully 
reproduces the overall features observed at the Tevatron. In this formulation, 
the average color content of the hard scattering is well-defined. We find that 
hard dijet rapidity gaps in quark-antiquark scattering are not due to singlet 
exchange alone. 
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Among the most intriguing recent experimental results in quantum chromodynam- 
ics is the observation of dijet rapidity-gap events, with anomalously low radiation in 
a wide interjet rapidity region [|. These events are typically identified by low or 
zero hadron multiplicity in the central region, despite the high momentum transfer 
necessary to produce the jets. 

The existence of such events was originally suggested on the basis of color flow 
considerations in QCD f|, ||. H forward jets are produced by exchanging a pair 
of gluons in a color singlet state, color can be recombined independently in each 
forward region. Then much less radiation is expected between the jets than when the 
exchange is a color octet gluon, which requires recombining color between particles 
moving in nearly opposite directions. Rapidity gap events have special interest as 
clear illustrations of color coherence and its interplay with hadronization. In addition, 
because their observation requires large rapidity intervals, they offer a new window 
into a perturbative, yet Regge-like limit of QCD. Nonetheless, despite their intuitive 
appeal, the theoretical understanding of rapidity gaps has been somewhat hampered 
by two problems. One of these is the issue of "survival" |5|, [|. In any high-energy 
scattering, multiple soft interactions between spectators of the hard interaction may 
fill the gap by processes unrelated to the color content of the hard interaction. The 
second is that, since even the softest gluon carries color in the octet representation, 
it is not immediately obvious how the color of the hard scattering is to be defined. 

In this paper, we observe that it is possible to overcome these problems, at least 
in part, by identifying rapidity gaps in terms of energy flow, rather than multiplicity. 
The energy flow Q c into the central rapidity interval between a pair of jets is an 
infrared safe observable. That is, da/dQ c can be written as a convolution of parton 
distributions with a perturbative hard-scattering function, which depends on Q c . The 
issue of color flow may then be formulated self-consistently in the hard-scattering 
function. Corrections to the factorized cross section are proportional to powers of 
A/Qc, with A the scale of the QCD coupling, and may become large for small Q c . 
As we shall see, however, the purely perturbative cross section remains well-defined, 
and energy flow gaps appear in this limit, once soft radiation is resummed including 
color effects [0, §, ||. In the conventional formulation for rapidity gaps, one writes 
/gap = /singlet -Ps, with / gap the fraction of gap events, / S i ng iet the fraction of "hard 
singlet" exchanges, and Ps the survival probability. Compared to this, we generalize 
/singiet, which then necessarily incorporates a perturbative survival probability. Non- 
perturbative survival considerations may reappear as we approach zero energy in the 
gap, but their importance should be reduced in a calorimetric measurement. Our 
results below support this possibility. 

To be specific, we will study the process p(pa) + P(pb) —> Ji(Pi) + Ji{P2) + A gap , 
where we sum inclusively over final states, while measuring the energy that flows into 
the intermediate region between two forward jets. For simplicity, we restrict ourselves 
to valence quarks and antiquarks, q{k,A) + <?(&b) — ► q{ki) + qfa) + X. We will begin 
by deriving a cross section for this process, specific to the geometry described by the 
DO and CDF collaborations |^, |TD|. We go on to evaluate the cross section as a 
function of Q c , and to view the results in the light of what we have learned from 
experiment. We will close with a few comments on the relation of our approach to 
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Figure 1: Geometry of the calorimeter detector. Q c is the energy flow into the central 
rapidity interval 



previous work, and on prospects for further development in this problem. 

Following CDF and DO, we require the two jets, and therefore the outgoing partons 
coming from the hard scattering, q(kx), 5(^2), to be directed into fixed forward and 
backward (collectively denoted "forward") regions of the calorimeter, defined by \y\ > 
yo, where y is the (pseudo)rapidity y = (1/2) lncot(6 l /2), with 9 the polar angle. 
In addition we require the jets to have transverse energies above an experimental 
threshold, Et- We will discuss cross sections for measured energy in a symmetric 
central region, spanning rapidity Ay = 2y$. This geometry is presented schematically 
in Fig. |l[ The inclusive dijet cross section for all events with energy in the central 
region equal to Q c is a typical factorizable jet cross section, which may be written as 

da 



dQ c 



(S,E T ,Ay) = ]T J d cos 6 

IaJb=u,cL 

x J dxA J Q dx B (f>f A /p( x A,-i) <frf B /p(xB,-i) 

^ da^ /~ A 

x 1^ -77—, j-[t,s,yjj,Ay,a 8 {t) 

f u f 2 =u,ddQ c d cos 9 



where <f>f A / p , 4>j B /p are valence parton distributions, evaluated at scale —t, the dijet 
momentum transfer, da^ / dQ c d cos 9 is a hard scattering function, starting with the 
Born cross section at lowest order. The index f denotes /a + /s — ► fi + fz- The 
detector geometry determines the phase space for the dijet total rapidity, yjj, the 
partonic center-of-mass (cm.) energy squared, s, and the partonic cm. scattering 
angle 6, with — | — cos^ = i. For simplicity of presentation, we take Q c to be the 
energy in the dijet cm. . 

In the spirit of Refs. 0, || |9] , we now observe that we may perform a further fac- 
torization on the partonic hard-scattering function da^/dQ c d cos 9. The underlying 

observation is that for Q c \T— the soft gluon radiation that appears in the central 
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region decouples from the dynamics of the hard interaction that produces the dijet 
event. In technical terms, soft gluon emission may be approximated by an effective 
cross section, in which the hard scattering is replaced by a product of recoilless color 
sources (specifically, Wilson lines || ||) in the directions of the incoming partons and 



the outgoing jets. The refactorized hard-scattering function then takes the form [11 



dQc v ' \ n 

xS LI (^,yjj,AyY (2) 

The functions Sli and Hjl contain the dynamics of soft radiation from Wilson lines 
(at measured Q c ), and the hard interaction, respectively, with \x a new factorization 
scale. The product itself must be independent of how we choose the scale, so long 
as fi > Q c . The indices I and L label the possible color structures of the hard 
interaction, which correspond in Sli to the color matrices that couple the four Wilson 
lines representing the 2 — > 2 hard subprocess. One index refers to the hard scattering 
in the amplitude, the other to the hard scattering in the complex conjugate. For 
quark- ant iquark scattering, the Wilson lines are in the 3 (quark) and 3* (antiquark) 
representations of SU(3), respectively, and their product may be characterized by 
either singlet or octet color exchange in the t- or s-channel. For the physical reasons 
outlined above, we will choose a t-channel color basis. Corrections to Eq. @ are 
expected from three-jet final states, for which the analysis below must, in principle, 
be repeated. We expect these corrections to be relatively small. 

Because the left-hand side of Eq. (0) is independent of the precise choice of factor- 
ization scale fi, the matrices Hjl and Sli must satisfy evolution equations, in which 
their variations with \x cancel each other. The only variable that H and S hold in 
common is a s (fi 2 ), and, as a result, the evolution equation for S is 

+ (3{9) ^) Su = -( T ^ lbSbi ~ SlaFsUi , (3) 

and similarly for H, with r,s(a s ) an anomalous dimension matrix. Consider a t- 
channel singlet-octet basis, with color vertices schematically given by 

Cl = /x J, c 2 = Y,T a xT a , (4) 

a 

with / the identity and T a the generators of SU(3) in the quark representation. A 



one-loop calculation in this basis gives [I I 



where iV c is the number of colors. The functions £ and p are 

N 2 — 2 

£(A|/) = _2iV c Ay + 2m^^, (6) 
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p(yjj, Ay, . 
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'cos(0) - tanh 



cos(#) + tanh (- 
2 
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2 A iV„ 2 
+— Ay - 2tn ' 



(7) 



While p depends on the jet rapidities and on the partonic scattering angle, £ depends 
on the geometry only, through Ay. We note that the off-diagonal components of Tg 
are purely imaginary. This interesting feature is due to strong coherence effects in 



the one- loop calculation, related to angular ordering QT3 1 . 

To study the Q c -dependence of S, it is convenient to diagonalize F$- In the basis in 
which Ts is diagonal, Eq. (^) implies that the components of S evolve independently 
in [i. In this basis we may calculate unambiguously the dependence on the central 
energy flow Q c . This is the technique that we summarize in the following. 

The eigenvectors of Ts in Eq. (|5|) may be chosen as 



ei 



( :, 2 



8^ + 7K 



-i 



where we define 

V (Ay) = V / iV c [£(A ? /)] 2 -32C F 7r 2 . (9) 

A very useful feature is that these eigenvectors are independent of the jet rapidities, 
and depend only on Ay. The corresponding eigenvalues of Ts are in general complex, 



Ai 



2tt 



1 



— p 

2 H 2, 



A 2 = ^ 



2~n 



1 1 
2 P+ 2- 



N, 



■V 



(10) 



In the limit of a large central region, Ay ^> 1, the function £ has a large negative 
real part, while the real part of r] is positive, and grows with Ay. From Eq. (|8|), 
we see that, as Ay — > oo, e\ reduces to a color "quasi-singlet" , and e 2 to a color 
"quasi-octet" . As a realistic example, we take the value Ay = 4, and find 



(=2 



0.455 e 2 - wli j 
0.101 e- 0981i 



(11) 
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For this configuration, the second eigenvector is close to a color octet, but the first 
is still a mixture of octet and singlet, with the latter only slightly predominant. 
In the following, however, we find it suggestive to retain the names "quasi-singlet" 
and "quasi-octet" for the elements of the diagonal basis. In the limit of large Ay, the 
eigenvalue for the quasi-octet grows with Ay, while the eigenvalue of the quasi-singlet 
does not. This will produce the expected enhancement of the latter relative to the 
former in the resummed cross section at small Q c . We shall use Greek indices to 
identify the basis in which r$ is diagonal. 

We can now write down a resummed cross section, working to lowest order in 
a s (—t), but resumming all leading logarithms in Q c . We transform Eq. @ to the 
diagonal basis, and solve the evolution equation for S, to get 



da^ 



dQ c d cos 9 



s,i, yjj,Ay,a s (-i) 



(Ay, V§, \f-t, a s (-t) ) fig (Ay) 



E 



x- 



7/9 



Qc 



In 



Qc 

A 



-"tp- 



ln 




—E, 
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(12) 



The coefficients E y p are given by 



2tt 



E lP (yjj, 6, Ay) = — [\; (yjj, 6, Ay) + \ p (yjj, 6, Ay 



(13) 



where (3\ is the first coefficient in the expansion of the QCD /^-function, f3\ = yA^ ( 
|n/, and where we define by = a, s \p + • • -. 



In accordance with our approximation, the matrix S^l is obtained by transforming 



the zeroth order 
traces, 



li of E q- 



to the new basis. The matrix S L °j is just a set of color 



o(o) 




\ (N? - 1) 



(14) 



and is transformed to the diagonal basis by the matrix (R 



>Kf3 



o(0) 



Analogously, we take for Hj2 the square of the single-gluon exchange amplitude, rep- 
resented in the color basis. Considering the dominant t-channel Born-level amplitude 
alone, which is purely octet, we have H IL = S12 <5l2 &t, where a t is the t-channel 
partonic cross section, including the coupling a s (—t). The contribution of s-channel 
diagrams has a relatively small effect, and will be described elsewhere In the 

diagonal basis the hard matrix Hjl becomes Hj$, defined as 



R- 



7 m 



O(0) I D- 



N/3 



(15) 



H 



(i) 



( R )(3L H LK 



(rf) 



(16) 
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E=630 GeV 



Qe/(1 GeV) 



E=1800 GeV 



Q c /(1 GeV) 



Figure 2: The cross section (solid line) and the contributions from quasi-octet (dotted 
line) and quasi-singlet (dashed line), for y/S = 630 GeV, Ay = 3.2, and \/S = 1800 GeV, 
Ay = 4.0, respectively. Compare Fig. 1 of Ref. [1C]. Units are arbitrary. 



Observe that and both acquire a Ay-dependence through the change of 
basis. 

From Eqs. fll2|)-(|13D, using the results described above, it is possible to evaluate Eq. 
(H]). For the valence partons we have taken the leading order CTEQ4L distributions 
12]. In Fig. ^] we plot the shapes of the cross sections obtained in this way, as a 



function of the radiation into the central region, Q c , for two different sets of conditions, 
= 630 GeV, Ay = 3.2, and = 1800 GeV, Ay = 4.0. We also show the 
contributions of quasi-octet and quasi-singlet terms. There is in addition a negative 
interference term, not exhibited separately. As anticipated above, we find a strong 
suppression of the quasi-octet component for very small values of Q c , contrasted to 
a peak for the quasi-singlet in this limit. The reason for this difference is easily 
found. For most kinematic configurations, the coefficient En from Eq. (f[3|) is less 
than one, so that the quasi-singlet cross section in Eq. (|l^) decreases monotonically 
with increasing Q c . For the quasi-octet, on the other hand, E 2 2 is always greater 
than unity, so that its contribution grows with Q c , until the power of the logarithm 
is overcome by the dimensional factor of 1/Q C - 

These results can be compared with the experimental data in Fig. 1 of Ref. ]10|, 



showing the measured number of events as a function of the number of towers counted 
in the central region of the calorimeter, clearly related to Q c . We can understand the 
similarity of shapes in terms of the Q c dependence in Eq. (|I2"D, discussed above. This 
similarity is suggestive; indeed, from our simulation we have evaluated the minimum- 
maximum ratio of the cross section, finding about 30% at y/S = 630 GeV and about 
15% at y/S = 1800 GeV, close to the analogous ratios in Fig. 1 of Ref. [|1(J. We have 
also determined an analog of a "hard singlet fraction" 0,11, as the ratio of the area 
under the quasi-singlet curve to the area under the overall curve. It is about 5% at 
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= 630 GeV and about 3% at VS = 1800 GeV. The order of magnitude of the 
result is reasonable, although higher than the roughly 1% found at the Tevatron using 
track or tower multiplicities. How much of this difference is due to our new defini- 
tion of the gap and how much to the lack of a nonperturbative survival probability 
remains to be explored. The sharp upturn that we observe below 1 GeV is due to the 
divergence of the perturbative running coupling at Q c = A; nonperturbative effects 
will attenuate this rise. 

Previous analysis of rapidity gaps in dijet events has tended to emphasize either 
the short- distance ^ or long-distance fl5], [16|, [L7| aspects of the problem. (The 
role of Sudakov logarithms in double-rapidity gap events has been discussed in |Tj§ .) 
Here, we have argued that by factorizing short- and long-distance effects, we may treat 
both dependences systematically. In our formalism, the mixing of color states begins 
at short distances precisely with two-gluon exchange || (L4[], summarized through the 
anomalous dimension Ts, while long-distance color ("bleaching") effects [T5|, |TB|, |TT| 
follow the evolution of the different color components between the short-distance scale 

and the long-distance scale Q c . 
As we have observed above, our formalism does not include a nonperturbative 
survival probability |5|, |6| associated with the interaction of spectator partons. Clearly, 
a full phenomenological analysis will also require the inclusion of processes involving 
gluons (including qq — > gg) and sea quarks. The treatment of gluon-gluon scattering 
fi~6| , |i~7f should be particularly interesting |J. Nevertheless, we believe that the basic 
features shown in the valence-quark analysis outlined above will appear as well in 
a more complete discussion. A calorimetric analysis of dijet rapidity gap events, if 
possible experimentally, could shed valuable light on the dynamics of QCD. 
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